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Abstract. Distortions of CMB temperature and polarization anisotropy maps caused 
by gravitational lensing, observable with high angular resolution and sensitivity, can 
be used to constrain the sterile neutrino mass, offering several advantages against the 
analysis based on the combination of CMB, LSS and Lya forest power spectra. As 
the gravitational lensing effect depends on the matter distribution, no assumption on 
light-to-mass bias is required. In addition, unlike the galaxy clustering and Lya forest 
power spectra, the projected gravitational potential power spectrum probes a larger 
range of angular scales, the non-linear corrections being required only at very small 
scales. 

Taking into account the changes in the time-temperature relation of the primordial 
plasma and the modification of the neutrino thermal potential, we compute the 
projected gravitational potential power spectrum and its correlation with the 
temperature in the presence of DM sterile neutrino. 

We show that the cosmological parameters are generally not biased when DM sterile 
neutrino is included. From this analysis we found a lower limit on DM sterile neutrino 
mass m„ s > 2.08 keV at 95% CL, consistent with the lower mass limit obtained from 
the combined analysis of CMB, SDSS 3D power spectrum and SDSS Lya forest power 
spectrum (m„ s > 1.7 keV). 

We conclude that although the information that can be obtained from lensing 
extraction is rather limited due to the high level of the lensing noise of Planck 
experiment, weak lensing of CMB offers a valuable alternative to constrain the dark 
matter sterile neutrino mass. 
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Introduction 

The latest results of the cosmic microwave background (CMB) temperature and 
polarization anisotropies from WMAP 3-year observations [U |2] combined with other 
tracers such as large-scale structure (LSS) galaxy surveys, supernovae luminosity 
distance and Lya forest, have lead in specifying the ACDM model as the cosmo logical 
concordance model [3j 0] . 

The direct confirmation of the ACDM theory is the detection of the acoustic Doppler 
peaks structure of the CMB angular power spectrum. Further successes are related to 
the correct prediction of the hierarchical structure formation via gravitational instability, 
the abundance of clusters at small redshifts, the spatial distribution and the number 
density of galaxies, the LSS matter power spectrum, the Lya forest amplitude and 
spectrum [5]. Despite its successes on large scales, the ACDM model produces too 
much power on small scales. In general, the observed structures have softer cores, lower 
concentrations and are less clumped than those predicted by the ACDM model (see Ref. 
[6] and references therein). 

A possibility to alleviate the accumulating contradiction between the ACDM model 
predictions on small scales and the observations is to add properties to the dark matter 
(DM) sector, relaxing the hypothesis on dark matter as being cold. Free streaming due 
to the thermal motion of the DM particles is the simplest known mechanism for smearing 
out small scale structure. Between the large free streaming distance of Hot Dark Matter 
(HDM) particles and the small free streaming distance of Cold Dark Matter (CDM) 
particles lies the intermediate scale of Warm Dark Matter (WDM) particles [7J [8]. 
The velocity dispersion of WDM particles is sufficient to alleviate some contradictions 
between the ACDM model predictions and observations such as: the predicted number of 
halos compared with the observed number of satellite galaxies in the Local Group [H E], 
the prediction of l/r a (1 < a < 1.2) behavior for galaxy rotation curves compared with 
the observed linearly rising behavior [TO"! [TTl [12], the prediction of too much baryonic 
material with low angular momentum to form the observed rotationally disk galaxies 

On the other hand, in addition to the evidence for mixing of active neutrinos 
from solar and atmospheric oscillation experiments, there are indications for another 
oscillation with larger mass-squared difference coming from short base-line neutrino 
oscillation experiments [TU [15] that can be explained by adding one or two sterile 
neutrinos with eV-scale mass to the standard scheme with three active neutrinos (see 
Ref. [16] for a recent analysis) . 

Such results have impact on cosmology because sterile neutrinos can contribute to the 
number of relativistic degrees of freedom at the Big Bang nucleosynthesis [17J. These 
models are subject to strong bounds on the sum of active neutrino masses from the 
combination of various cosmological data sets [TU [20] that rule out a thermalized sterile 
neutrino component with eV mass-scale [41 IT9]. 

However, there is the possibility to accommodate the cosmological observations with 



Constraints on non-thermal Dark Matter from Planck lensing extraction 



3 



data from short base-line neutrino oscillation experiments by postulating that sterile 
neutrino is not thermalized and has a phase-space distribution significantly suppressed 
relative to the thermal distribution. It is already known that sterile neutrino with 
mass of few keV - the Tremaine-Gunn [21] bound - provides a valuable DM candidate 
[221 1231 EH E5]. In the standard non- resonant production mechanism (small lepton 
number of the order of baryon number, L ~ 10~ 10 ) sterile neutrinos are produced via 
small mixing angle oscillation conversion of thermal active neutrinos [26j . 

Sterile neutrino with the mass in keV range has a radiative decay channel emitting 
a photon with an energy that is half of its mass eigenstate. The width of the decay line 
increases as the fifth power of its mass eigenstate and as square of its mixing angle, being 
potentially detectable in various X-ray spectra of the astrophysical objects [271 [28]. The 
strategy to search for dark matter particles possesing a radiative decay channel and to 
derive the constraints on their parameters from X-ray observations applied to sterile 
neutrino is discussed in Ref . [291 EO] • 

The analysis of the observed X-ray background from HEAO-1 and XMM-Newton [31] 
leads to an upper mass limit for sterile neutrino of m Vs < 8.9 keV [jj improved to 
m Vs < 6.3 keV from the combined analysis of XMM-Newton observations of the Virgo 
and Coma clusters [32]. Also, the constraints on the rate of sterile neutrino radiative 
decay obtained from the analysis of the diffuse X-ray spectrum of Andromeda galaxy 
leads to m Us < 3.5 keV, which is a significant improvement over previous upper limits 
[33]. 

Although prone to uncertainties due to the estimate of DM distribution in dwarf galaxy 
[34] , interesting constraints on the parameters of radiative DM decay are obtained from 
XMM-Newton observations of the DM halo of the Milky Way and Ursa Minor [36]. The 
recent spectral analysis of the unresolved component of the cosmic X-Ray background 
in the Chandra North and South Deep Fields provides limits on the sterile neutrino 
mass [37] . The highest Milky Way halo mass estimate provides a limit on sterile neutrino 
mass of m Us < 2.9 keV in Dodelson-Widrow production model [22], while the lowest halo 
mass estimates provides a more conservative limit of m Us < 5.7 keV. 

The radiative decay of the sterile neutrino can also augment the ionization fraction 
of the primordial gas at high redshifts leading to the increase of the temperature of 
the primordial gas, the enhancement molecular hydrogen formation and of the star 
formation rate [381 GSJ H3 • 

The combined analysis of the Lya forest power spectrum measured by SDSS, the 
CMB anisotropy and the galaxy clustering power spectra yielded to the lower limits 
for sterile neutrino mass, m Vs > 14 keV [H] and m Vs > 10 keV [42] in Dodelson- 
Widrow production model [22], excluding sterile neutrino as DM candidate. However, 
the method based on the combined analysis of angular power spectra gives direct limits 
for the free streaming lengths of DM particles, the limits on their masses depending on 
their momentum distribution functions and therefore on their production mechanisms. 

|| Throughout the paper the sterile neutrino mass is quoted at 95% CL. 
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Taking into account the deviations from the thermal spectrum of DM particles 
produced due to the variation of the number of degrees of freedom leading to changes 
in the time-temperature relation of the primordial plasma and the modification of the 
neutrino thermal potential [2U [261 H3], the combination of CMB, LSS and Lya forest 
angular power spectra leads to m Vs > 1.7 keV with a further improvement to m Us > 3 
keV when high-resolution Lya forest power spectra are considered [4"3~] . 

In this paper we explore an alternative method to constrain the sterile neutrino 
mass based on the CMB weak lensing extraction. Distortions of CMB temperature and 
polarization maps caused by the gravitational lensing potential, observable with high 
angular resolution and sensitivity, have impact on cosmological parameter degeneracies 
when non-minimal cosmological scenarios are considered [4*4"l |4"5| l4*6l W7\ . 
The CMB weak lensing offers several advantages against the method based on the 
combination of angular power spectra. As the gravitational lensing effect depends on 
the dark matter distribution in the Universe, no assumption on light-to-mass bias is 
required. The projected gravitational potential is sensitive to the matter distribution 
out to high redshifts, preventing from non-linear corrections required only at very small 
scales. In addition, unlike the galaxy clustering and the Lya forest, the projected 
gravitational potential probes a larger range of angular scales, most of the signature 
comming from large scales. 

The paper is organized as follows. In Section 1 we compute the energy density 
evolution of the dark matter sterile neutrino in the expanding Universe, incorporating 
the time-temperature relation and the modification of the neutrino thermal potential. 
In Section 2 we compute the deflection angle power spectrum and its cross-corelation 
with the temperature when sterile neutrino dark matter energy density is considered. 
In Section 3 we derive limits on the sterile neutrino mass from Planck weak lensing 
extraction. We draw our main conclusions in Section 4. 

Throughout we assume a background cosmology consistent with the most recent 
cosmological measurements [3] with energy density of Q m = 0.3 in matter, = 0.05 in 
baryons, Q u = 0.01 in three active neutrino flavors, Qa = 0.7 in cosmological constant, 
a Hubble constant of H =72 km s -1 Mpc _1 , a reionization redshift z re =13 (assuming 
a sharp reionization history) and adiabatic initial conditions with a power-law scalar 
spectral index n s = 1. 

1. Sterile neutrino dark matter production and density evolution 

1.1. Non-resonant production 

Sterile neutrinos dark matter candidates are produced in the early Universe through 
oscillation with active neutrinos, at temperatures close to QCD phase transition [24, 26J. 
The oscillation process takes place due to the fact that the neutrino mass eigenstate 
components propagate differently as they have different energies, momenta and masses. 
Eigenstates of neutrino interaction include the active neutrinos, v a (a=e, /i, r), which 
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are created and destroyed in the standard model by weak interactions as well as the 
sterile neutrinos, u s , which do not participate in weak interactions. For oscillations 
occurring in vacuum between two neutrino flavors u a and v Sl the mixing can be written 
as: 

\v a >= cos 6*o |^i > +sin6'o|z/2 > 

\v s >=— sin 8 \i>i >+ cos 9 \is 2 > , (1) 

where v\ and v 2 are the neutrino mass eigenstate components and 8q is the vacuum 
mixing angle. The mixing of antineutrinos can be obtained from the above equation 
by performing the transformation v a —>■ v a and v s — > v s . It is usual to consider that 
each neutrino/antineutrino of a definite flavor is dominantly one mass eigenstate. In 
this circumstance we refer to the dominant mass eigenstate component of v a jv a as vi 
to that of v s jv s as v 2 and to their difference of squared masses as d~m 2 = m 2 Vs — m 2 . 

The mixing of the mass eigenstate components are modified in the presence of 
the finite temperature background either in the case of small active neutrino lepton 
number, L Ua ~ 10 -10 (of the order of baryon number), driving non-rezonant sterile 
neutrino production [22], or in the case of several orders of magnitude larger lepton 
number [2U HE]. The lepton number is defined as the difference between the active 
neutrino and antineutrino number densities normalized by the photon number density: 
L Va = {n Va -n 5 J/n 7 . 

The matter mixing angle, 8m, is related to the vacuum mixing angle, 8 , through 
(see e.g. Ref. [24] and references therein): 

■2 2f) = A 2 (p)sm 2 26 

M A 2 (p)sm 2 29 + D 2 ( P ) + [A( P )cos29 + V L -V T ( P )} 2 ' {> 

In the above equation A(p) = 5m 2 /2p is the vacuum oscillation factor, V L is the 
asymmetric lepton potential, V T (p) is the thermal potential, D{j>) = T Ua (p)/2 is the 
quantum damping rate and T Va is the collision rate defined as: 

1.27G 2 F pT\ a = e 
Q.92G 2 F pT 4 , a = h,t 



The asymmetric lepton potential, V L , and the thermal potential V T (p) read as [2 
V L = V2G 



2{n Ua - n Pa ) + 



V T = "^^« E Va > n Ua+ < E Pa > n Pa ] 
8V2G F p 



3m 



2 



< E Va > n Va + < Ep a > n 9a ) , (3) 

"w 

where: E Va is the active neutrino total energy, n Ua /n Pa are neutrino/antineutrino number 
densities, cosmological n n is the baryon number density, mz is Z° boson mass and mw 
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Figure 1. Panel a): The evolution with the scale factor of the temperature obtained 
by the integration of the time-temperature relation considering the variation of the 
statistical weight in relativistic particles g* , compared with the linear evolution of the 
temperature, Tu n = Tca/b(1 + z). Panel b): The evolution with the temperature 
of the statistical weight in relativistic particles g* used in the computation of the 
time-tcmpcraturc relation. Panel c): The temperature dependence on sterile neutrino 
comoving momentum obtained for the same sterile neutrino mass values as in Panel a). 



is W ± boson mass. The number density of the thermally distributed active neutrinos, 
n Va is given by: 

- h ft)' • < 4 > 

where n 7 = 2£(3)T^ /n 2 is the photon number density and £(3) is the Riemann zeta 
function of 3. 

Throughout we will consider a lepton symmetric Universe (n„ a ~ n„ a ) and neglect the 
contribution to the asymmetric potential due to the baryon number which is very small 
at the temperatures of interest [26]. In this case the contribution of asymmetric lepton 
potential V L in eq. (2) is negligible. 

1.2. Energy density evolution 



Since sterile neutrinos/antineutrinos are produced non-thermally, their mean energy 
density and pressure must be computed by the direct integration of their phase-space 
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Figure 2. Left: The dependence on the scale factor of the energy density parameters 
of sterile neutrinos expressed in units of the mean energy density of one mass eigenstate 
of massless neutrinos for different sterile neutrino mass values. Right: The evolution 
with the scale factor of the energy density parameters for matter, radiation and vacuum 
obtained for fiducial model when the corresponding sterile neutrino energy density 
parameters are included. 



distributions: 

rp4 poo 

pv s + pp 3 = / dq q 2 E Vs [f Vs (q) + f Ps (q)] , 
T 4 f 00 q 2 

where: f Ug and f Pg are the neutrino/antineutrino phase-space distribution functions and 
E Vs = A y 1 q 2 + a 2 m1 s is the total energy of sterile neutrino having a mass m Us and a 
comoving momentum q = E/T; a = (1 + z)^ 1 is the scale factor (ao=l today). The 
temporal evolution of sterile neutrino energy density affects the Hubble expansion rate 
that reads as: 

H{a) = ^^[Q m /a 3 + Q r /a 4 + Q A f 2 . (6) 

In the above equation G is the gravitational constant, Q m = Qb + ^cdm + &v a + is the 
matter energy density parameter where Q&, fl c dm, Q» a and Q Ug are the energy density 
parameters for baryons, cold dark matter, active and sterile neutrinos respectively, Q r is 
the radiation energy density parameter and fl\ is the vacuum energy density parameter. 



Constraints on non-thermal Dark Matter from Planck lensing extraction 



8 



Depending on the mass, the maximum rate production of sterile neutrino occurs at 
H3J: 

T_«133MeV(^) . 

For sterile neutrino with mass in the range of interest, the time evolution of the 
temperature-dependent thermal potential and of the collision rate require the knowledge 
of the time-temperature relation in the expanding Universe: 



da dptot 



(ptot+Ptot) \ (7) 



dT dT 

where p to t and p tot are the total density and pressure, including the contribution of 
sterile neutrino/antineutrino as given in equation (5). The time-temperature relation 
takes into account the variation of statistical weight in relativistic particles that is 
changed by nearly an order of magnitude since T max [2U H3]. A general treatment of 
the time-temperature relation, whose approach we have incorporated, is given in the 
Appendix of Ref. [24] . 

Left panel of Figure 1 presents the evolution with the scale factor of the temperature 
obtained by the intergration of the time-temperature relation taking into account the 
variation of the statistical weight in relativistic particles g*, compared with the linear 
evolution of the temperature Tu n = TcmbO- + z). We also show, in the right panel 
of the same figure, the dependence of the temperature on sterile neutrino comoving 
momentum. 

The time evolution of the sterile neutrino phase-space distribution, f Us , can be 
described by the semiclasical Boltzmann equation [221 [24] : 



j t UM) - Hit) q « i> a ^ v s )\UM - . ( 8 ) 

where: dt = da/aH(a), f Ua is the active neutrino phase-space distribution and T(v a ^ 
v s ) is the effective production/annihilation rate of sterile neutrinos: 

r(z/« ^ v s ) « Q.25T Va (p,T) sm 2 29 M . (9) 

We use the set of equations (5)-(9) to compute the temporal evolution of the energy 
density of sterile neutrinos in the expanding Universe. For sterile neutrino with the 
mass in the range of interest, the flavor dependence is almost negligible [26], implying 
that similar results can be obtained for all v a ^ v s mixings. Throughout we consider 
v T ^ v s mixing. 

Left panel of Figure 2 presents the dependence on the scale factor of the energy density 
parameter of sterile neutrino expressed in units of the mean energy density of one 
mass eigenstate of massless neutrinos for three different mass values. We also show the 
same dependence for the case m Vs = 0. The right panel of the same figure shows the 
evolution with the scale factor of the energy density parameters for matter, radiation 
and vacuum obtained for our flat ACDM background cosmology when the corresponding 
sterile neutrino masses are included. The production of DM sterile neutrino affects the 
radiation and matter evolution through the Hubble expansion rate, altering the time- 
temperature relation and redshifting the relativistic species. 
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2. CMB lensing extraction with Planck experiment 
2.1. Deflection angle power spectrum 

The CMB photons from the last scattering surface are subject to the cumulative effects 
of the large scale structure gravitational potential [50] • The net result is that the CMB 
temperature anisotropy and polarization patterns measured in the direction n in the sky 
provides information on the photons emerging from the last scattering surface from the 
direction n = n + Sn. The intensity and the linear polarization of the lensed CMB are 
completely specified by the Stokes parameters /, Q and U which are related with the 
unlensed Stokes parameters through: X(n) = X(h), where X(n) stands for the lensed 
/, Q and U. The deflection field is defined as the gradient of the lensing potential 
5n = V\l/(n), where V is the covariant derivative on the sphere. 

The deflection map and the lesnsing potential map can be expanded in spherical 
harmonics [511 [52], so that the relation between the deflection angle power spectrum 
Cf d and the lensing potential power spectrum Cf is given by [53J: 

Cf = 1(1 + l)Cf * . (10) 

The power spectra of the lensing potential, Cf and the correlation with the 
temperature anisotropy, Cf T , can be numeriacally computed by using the Code for 
Anisotropies in the Microwave Background (CAMB) [501 1521 151] in the linear theory 
as well as by including the non-linear corrections from HALOFIT [55]. We recall that 
Cf T does not vanish because the temperature map includes information on the time 
variation of the gravitational potential though the integrated Sachs-Wolfe effect. 
We modified the CMB anisotropy code CAMB to compute the CMB angular power 
spectra in the presence of a DM sterile neutrino component. We include in 
the computation the momentum-dependent sterile neutrino phase-space distribution 
function, as described in the previous section, its unperturbed and perturbed energy 
density and pressure, energy flux and shear stress. 

Figure 3 presents in the top panel the deflection angle power spectrum, Cf d , 
obtained in our fiducial model imv s = 0). The bottom panel of Figure 3 shows the 
relative percentual differences between the fiducial model and exact the same model but 
with m v = 1.7, 3.5 and 8.3 keV. 

The signature of sterile neutrino mass on the deflection angle power spectrum is the 
suppression of its amplitude relative to the fiducial model. The net suppression is scale 
dependent and the relevant length scale is the free-streaming scale [56] , which increases 
with the mean dark matter sterile neutrino velocity and decreases with its mass. 

3. Constraints on sterile neutrino mass from PLANCK lensing extraction 

3.1. Gaussian Likelihood function 

In the standard inflationary cosmology, the spherical harmonic coefficients of the 
unlensed temperature and polarization field in the sky obey the the Gaussian statistics 
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Figure 3. Top panel: The deflection angle power spectrum, Cf d , obtained in the our 
fiducial model (m Vs — 0) and the minimum variance noise power spectrum (dashed 
line) for Planck experimental characteristics given in Table 1. Bottom panel: The 
relative percentual differences between the fiducial model (m Ug —0) and exactly the 
same model with m„ s = 1.7, 3.5 and 8.3 keV. 
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Figure 4. The CMB angular power spectra for the fiducial model. The error 
bars include the Planck experiment instrumental noise and the cosmic variance. For 
graphical representation convenience, the error bars assume a multipole binning of 
Al = 7 until I - 70 and then Al ~ Z/10 [SU [63] ■ 
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and can be used to define the estimators of the covariance: 



21 + 



X lm a,be{T,E,B} 



(11) 



where Wi m and Xi m are the harmonic coefficients of the CMB temperature and 
polarization anisotropy maps. The corresponding full-sky likelihood function is given 
by: 

X 2 = -21n£(C|C) = (2/+ l^T^CC" 1 ] +ln|C|}, (12) 

where C is the covariance matrix of the observations and C is the corresponding 
covariance matrix of estimators. The covariance matrix of the observations reads as: 



/ Cj T + Nf T 

\ 



C EE + N EE 








\ 



C BB + N BB 



(13) 



In the above equation Cj , a,b — {T, E, B} are the unlensed power spectra of primary 
anisotropies, Nf* are the corresponding instrumental noise variances and the parity 
invariance (C BT = C BE = 0) it is assumed. 

Weak lensing correlates the lensed multipoles and the lensed sky is not any more 
Gaussian [57]. However, as shown in Ref.jHH], at least up to Planck resolutions and 
sensitivities, equation (11) remains approximately correct if the B-mode, that is noise 
dominated, is not included in the covariance matrix. The non-Gaussian corrections 
can be however important for the future high sensitivity CMB polarization experiments 
having a signal dominated B-mode power spectrum [59J. 

In this paper we consider a fiducial model with no significant amplitude of the primordial 
gravitational waves, omitting the B-mode from the parameter estimation analysis. In 
this case, the data covariance matrix read as: 



\ 



(14) 



where Cf ', a, b 



TT + N TT QTE 

C TE C EE + N EE o 

y Cf d Cf A + N? d J 

{T, E}, are the lensed CMB power spectra, iV" 6 are the corresponding 
detector noise power spectra, Cf is deflection angle power spectrum, jVj the noise 
power spectrum associated to the lensing extraction and Cf d is the power spectrum of 
the cross-correlation between the temperature and deflection angle. 
The weak lensing does not change the total variance in the CMB temperature anisotropy 
and polarization maps and hence lensed observations have the same variance as if there 
were no lensing [58, 59J. 

For each frequency channel we consider an homogeneous detector noise with the 
power spectrum, iVf^ given by [60] : 



N, 



ab 



,A c ) 2 exp^ + W ln2 , 



(15) 
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Table 1. The expected experimental characteristics for the Planck frequency 
channels considered in this work [61j : v is the frequency of the channel, 6^ is the 
FWHM, At and Ap are the sensitivities per pixel for temperature and polarization 
maps. 



V 


FWHM 




A P 


(GHz) 


(arc-minutes) 


(a* k) 


(A* k) 


100 


9.5 


6.8 


10.9 


143 


7.1 


6.0 


11.4 


217 


5.0 


13.1 


26.7 



where v is the frequency of the channel, 0& is the FWHM of the beam and A c , where 
c G (T, P), are the sensitivities per pixel for the temperature and polarization maps. 
The global noise of the experiment, N* , can be written as: 

-i 



N, 



ab 



E 

V 



ab\-l 



(16) 



The expected experimental performances of the Planck frequency channels considered 
in this paper [61] are presented in Table 1. 

The lensing noise power spectrum, Nf d , was numerically computed by using the 
minimum variance quadratic estimator of Okamoto and Hu [62]. By definition, the 
quadratic estimator is build from pairs (a, b) of observed temperature or polarization 
modes and its multipoles are given in the quadratic form: 

,M 



d(a, b) 



M 
L 



N 



;i7) 



IV mm! 



where a™ and 6™ correspond to the lensed CMB fields. The minimum variance estimator 
is obtained by finding the weights @i \ {L) that minimize the Gaussian variance: 



< d(a, b)fd(a\ [Cf + K ba b 



For the minimum variance estimator the lensing noise reads as: 

-i 



N, 



dd 



aba b \— 1 



aba b 



(18) 



(19) 



In the top panel of Figure 3 we show the lensing noise power spectrum, Nf d , for Planck 

experiment, obtained by using the minimum variance estimator method. 

Figure 4 presents the CMB angular power spectra for the fiducial model. The error 



bars, AC", include the instrumental noise and the cosmic variance: 



AC? 



(21 + l)Alf sky 



- (C? + a= [TT, EE, TE, dd, Td] . 



(20) 
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Figure 5. Left: The marginalized likelihood probabilities for Planck- like 
observations obtained by using the lensing extraction with (solid red line) and without 
(dashed black line) DM sterile neutrino component. Right: The marginalized likelihood 
probabilities for the (inverse) mass of DM sterile neutrino obtained for the Planck- like 
observations with lensing extraction (solid red line) compared with the marginalized 
likelihoods obtained from [43]: CMB plus SDSS 3D power spectrum (dashed line), plus 
SDSS Lya forest power spectrum (dotted line), plus VHS high-resolution Lya forest 
power spectrum (dash-dotted line). 



In the above equation f s k y represents the fraction of the sky covered by the observations. 
For the purpose of this work we take f s ky = 0-S and assume a perfect cleaning of all 
the astrophysical foregrounds. For graphical representation convenience, the error bars 
presented in Figure 4 assume a multipole binning of Al = 7 until I ~ 70 and then 
Al ~ Z/10 [HHE3]: 

3.2. Parameters estimation 

We assume a flat ACDM cosmological model with the following set of parameters to be 
determined from the PLANCK-like observations: 

© = (Q b h 2 ,Q c h 2 ,9 s ,z re , f u ,m Ua ,n s ) 

where Qf, and fl c are the fractions of critical density in baryons and cold dark matter, 
9 S is the angular acoustic peak scale of the CMB, z re is the reionization redshift (we 
assume a sharp reionization history), /„ is the fraction of massive neutrinos, m Vs is the 
sterile neutrino mass and n s is the spectral index of primordial adiabatic perturbations. 

We use the CosmoMC [M] Monte Carlo Markov Chain (MCMC) public package 
with the modification of the function ChiSqExact [63] required to compute the 
likelihood function when the deflection angle power spectrum is included, to sample 
for the posterior distribution of parameters space 0, giving the PLANCK-like data. For 
graphical representation convenience, the error bars in Figure 4 assume a multipole 
binning of Al = 7 until I ~ 70 and then Al ~ Z/10 [fill 1KB]: 
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Left panel from Figure 5 presents the marginalized likelihood probabilities of the 
cosmological parameters for PLANCK-like observations obtained from lensing extraction 
with and without the DM sterile neutrino component, showing that the cosmological 
parameters are generally not biased when including DM sterile neutrino. 
In the right panel of Figure 5 we present the marginalized likelihood probability for 
the inverse mass, mj 1 , of sterile neutrino obtained for PLANCK-like observations by 
using the lensing extraction. We find a lower limit on the sterile neutrino dark matter 
m Vs > 2.08 keV at 95% CL, in agreement with our previous result [65J obtained from 
Planck lensing extraction by using the Fisher information matrix (m Us > 1.75 keV). 
For comparison we also show in the same figure the marginalized likelihoods for m" 1 
obtained in Ref. [13] from: CMB plus SDSS 3D power spectrum, plus SDSS Lya forest 
power spectrum, plus VHS high-resolution Lya forest power spectrum. Our result is in 
agreement with the result from the combined analysis of CMB, SDSS 3D power spectrum 
and SDSS Lya forest power spectrum (m Va > 1.7 keV) but is less constrained than the 
limit obtained when VHS high-resolution Lya forest data are included (m Us > 3 keV). 

4. Conclusions 

From the future high precision observations of the Planck satellite we would like 
to extract information about both unlensed and lensed CMB anisotropics. The weak 
lensing effects of the CMB induced by the neighbouring galaxy clustering will have a 
significant effect on the statistics of the observed CMB and must be take into account 
in order to obtain reliable parameter estimates, offering several advantages against the 
cosmological parameter extraction based on the combination of CMB, LSS and Lya 
power spectra. As the gravitational lensing effect depends on the matter distribution, no 
assumption on light-to-mass bias is required. In addition, unlike the galaxy clustering 
and Lya forest power spectra, the projected gravitational potential power spectrum 
probes a larger range of angular scales, the non-linear corrections being required only 
on very small scales. 

However, the amount of information that can be obtained from the CMB weak 
lensing extraction depends on the lensing noise level. 

In this paper we address the possibility to constrain the DM sterile neutrino mass from 
PLANCK-like observations by using the CMB lensing extraction. We use the minimum 
variance estimator of Okamoto and Hu [62] to compute the lensing noise power spectrum 
for Planck experimental configuration. 

Taking into account the changes in the time-temperature relation of the primordial 
plasma and the modification of the neutrino thermal potential, we compute the projected 
gravitational potential power spectrum and its correlation with the temperature 
in the presence of DM sterile neutrino componet that together with the lensed 
CMB temperature and polarization anisotropy power spectra are used to constrain 
the cosmological parameters including sterile neutrino mass, from PLANCK-like 
observations. 
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We show that the cosmological parameters are generally not biased when DM sterile 
neutrino is included. From this analysis we found a lower limit on DM sterile neutrino 
mass m Vs > 2.08 keV at 95% CL, consistent with the lower mass limit obtained from 
the combined analysis of CMB, SDSS 3D power spectrum and SDSS Lya forest power 
spectrum (m Vs > 1.7 keV) [43J. 

We conclude that, although the information that can be extracted is rather limited due 
to the high level of the lensing noise of Planck experiment, the weak lensing of CMB 
offers a valuable alternative to constrain the dark matter sterile neutrino mass. 
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